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14. Alihtct 

Close, to stationary periodic motions of an axisyrametjric satel- 
lite In a circular orbit are considered. The satellite was 
equipped with a spherical magnetic damper.' The Investigation 
was conducted on the assumption that a strong magnet was In- 
stalled on the damper float* Stationary rotations of the 
satellite around the axis of symmetry ’are* selected as the 
generating solutions; The solutions. are constructed In the 
form of .power series of the small parameter, and they are ex- 
tended numerically to the region of random values .of. the dajnp- 
ing coefficient. The stability of the resulting solutions was 
investigated. 
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ANNOTATION 


Close to stationary periodic motions of an axisymmetrlo satel- 
lite In a circular orbit are considered. The satellite was equipped 
with a spherical magnetic damper. The investigation was conducted 
on the assumption that a strong magnet was Installed on the damper 
float. Stationary rotations of the satellite around the axis of sym- 
metry are selected as the generating solutions. The solutions are 
constructed in the form of power series of the small parameter, and 
they are extended numerically to the region of random values of the 
damping coefficient. The stability of the resulting solutions was 
investigated. 
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PERIODIC MOTIONS (CLOSE TO STATIONARY) OF AN AXISYMMETRIC 
SATELLITE WITH MAGNETIC DAMPING 

M.Yu. Ovchinnikov 


1. Introduction 


Periodic motions of an axisymmetric satellite equipped with a /3* 
spherical magnetic damper, in orbits of random declination, were in- 
vestigated in [1]. With A>V7 (A is the ratio of the axial moment of 
inertia of the satellite to its equatorial moment of inertia), in- 
stability of motion of the satellite in the plane of the polar orbit 
was found with respect to spatial perturbations. An example of a 
transition process was presented there which resulted in stable mo- 
tion, which was characterized by deviation of the axis of symmetry of 
the satellite from the plane of the polar orbit by "small tremors" 
about this position and rotation of the satellite around the axis, 
with a period close to the period of rotation of its center of mass 
around the orbit. It is shown in the present study that, with a mag- 
netic damper aboard the satellite, its stationary rotations change to 
forced periodic motions. In the (sin 2 i, A) plane (i is the declina- 
tion of the orbit of the center of mass of the satellite to the plane 
of the equator), regions of existence of stationary rotations of the 
satellite are constructed, which are selected as generating motions. 
Motions of the satellite which are close to stationary rotations were 
constructed in the form of a power series of the small parameter. The 
ratio of the characteristic values of the damping and gravitational 
moments acting on the satellite was used as the small parameter. The 
orbit of the center of mass of the satellite is considered circular. 

The geomagnetic field is approximated by the field of a dipole which 
coincides with the axis of rotation of the earth. The resulting mo- 
tions can be used as nominal (operating) motions of an axisymmetric 

■Numbers in the margin indicate pagination in the foreign text. 
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satellite with X<1. 

*W 

The stationary rotational notions of an axlsynnetrlo satellite 
with model damping were constructed and studied In C2]. 

V.A. Sarychev and Yu. A. Sadov are thanked for attention to the 

work. 

2. Equations of Motion and Formulation of the Problem 

We consider that the satellite Is a solid body, the moments of 
inertia of the damper float are negligibly small compared with the 
moments of Inertia of the satellite, and the center of mass of the 
float Is fixed relative to the satellite. The motion of the float 
about Its center of mass then does not affect the Inertial character- 
istics of the satellite, and the float Is replaced by an equivalent /4 
point mass in their determination. 

To write the equations of motion of the satellite and float 
relative to the center of mass, we introduce the following clock- 
wise rectangular coordinate systems: 

OXjXgX^ bound to the satellite coordinate system; Its 
axes are the principal central axes of Inertia of the satellite; 
point 0 Is the center of mass of the satellite; 

OX^X^ 18 the orbital coordinate system; the OX^ axis is 
directed along the radius vector of point 0 relative to the center of 
mass of the earth; the OX^ axis coincides with the transversal, and 
the 0X 2 axis coincides with the normal to the plane of the orbit; 

Is the unit vector of the 0X 2 axis; 

0 Zl Z 2 Z 3 Is the magnetic coordinate system; the OZ^ axis Is 
directed along vector ft of the geomagnetic field strength at point 0; 
unit vectors 
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determine the corresponding axes of the magnetic coordinate system. 

We assign the position of the Ox^XgX^ coordinate system relative 
to the orbital coordinate system by means of angles a, 6» y (Fig. 1). 
The transition matrix and its elements have the form 



X* X A 

X/ 

ct„ °** 

X* 

<*33 

*9 

&jf Qti 


crtt*co*doosM , 

a* • -cosJL s/nB t 
Of* * $/**, 


Ofr*- s/nf coiaL, 

a M • sin f sin/3 -cctfsind oas/t, 
■s//y f ocs/$ toctft/tt* *’»/*, 
a M *c0sfcpsd. 


We assign the position of vector I of the magnetic moment of the 
magnet Installed in the float relative to the magnetic coordinate s s- 
tem by means of angles and The corresponding directing cosines 
have the form 

c f - coSol t cos p, y c^ss/njSf 9 Cj * -sinoi t ces . 


Let arbitrary vector q be assigned by projections q^qjtQj on 
the axes of any of the coordinate systems introduced , Ox^XgX^ for 
example. We will then write q*(q^»42» q 3^ x * etc * As nee<led » we will /5 
gmeley summation over the recurrent Indices and free indices. The 
indices run through the values 1, 2 and 3* 


We approximate the geomagnetic field by the field of a magnetic 
dipole placed at the center of the earth, the axis of which is directed 
along its axis of rotation. The projections of vector 3 at point 0 on 
the orbital coordinate system axes, referred to the quantity H v «um/p^, 
have the form 

Hf *V/tic0SU, Hgzccsc , 5//7 i 
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where um«8. 06*10^ Oe^cnr 1 1b the negnetlo moment ^of^e dipole, p it 
the radius of the satellite orbit, 1 is its declination to the plane 
of the equator, and u is the argument of the latitude. We assign the 
position of the magnetlo coordinate system relative to the orbital 
coordinate system by means of the transition matrix 


where 



2 / Z * h. 

hff h* hi$ 

h u h u 

hjj hjt 






Of the external moments which act on the satellite, we will take 
only the gravitational moment into account* Of the external moments 
which act on the damper float, we will take into account only the mag- 
netic moment. The interaction of the satellite and the float is due 
to eddy currents Induced in the outer shell of the damper by the mag- 
netic field of the float. The hypothesis was introduced above that 
the moments of inertia of the float are significantly less than the 
moments of inertia of the satellite. Therefore, with the exception of 
the small time Interval after freeing the float, its motion is de- 
termined by the dynamic equilibrium of the magnetic moment and the 
moment of the induced eddy currents [33. The motion of the axlsym- 
metric satellite and the damper float are described by the system of 
equations [1] 

: (i) 

& * - hj§c*)t 
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Here 

Jf* 

* %<&**•*>*•*>*'”<■ i - zu xu 

m,rnO, fiftttj,**— 

<-i&, «/•» 

A, C are the equatorial and axial momenta of Inertia of the satellite, 
wq la the angular velocity of orbital motion of the satellite center 
of mass, la the damping coefficient, W^W^n are the projections, 
related to of the absolute angular velocity of the satellite on 
the Resal* axes, which coincide with the axes of coordinate system 
Ox^XgX^ at 0*0. The Ox^ axis Is the axis of symmetry of the satel- 
lite. The point Is designated by differentiation over u. The deriva- 
tion of system of Eq. (1) Is described In [1], where another sequence 
of flight angles a, 0, y Is Introduced. 

Letl/e«I. The motion of the satellite Is then [1] described to 
within 0(I/e) by the following equations: 

a; r-M- HlfL* 

.jj. - . , cc$tL ] , ( 2 ) 

rh< f *nrM!(huC*fi 

Q * hj,ru}Jhj S ) +4u-&] , 


Here, 


A* # W f IhitCcU 

+ Ws(hfOOMr+h f Stof)+IlttM ¥ 0*, , 


Ve will next Investigate system (2). Angle 0, which describes the 
rotation of the satellite around Its axis of symmetry. Is determined 
by the equation 
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System of Eq« (2) with 2ir periodic clockwise segments with respect 
to u contains three parameters 1, A, kg. The latter two satisfy the 
Inequalities 0<A{2, kg>0. Equation- (2) do- not change their form In 
the following substitutions of the phase variables and parameter 1: 


a [iZ ?? ; 

w, — w,, 


(3) 


In accordance with this. It Is sufficient to Investigate the solu- 
tions of Eq. (2) In the Interval 0<i«V2. 

With kg*0, system (2) permits solutions which correspond to 
stationary rotations of the satellite. With kg>0, we obtained and 
Investigate the 2 ir periodic with respect to u solutions of system 
(2) generated from them. 

3. Stationary Rotations of a Satellite with kg«0 

If kg*0, system (2) has the generalized energy Integral 

(H) 

and stationary solutions y"Yo* a*«o OQ^const), which are de- 
termined by the system of equations 

cos ft [&&, +l i *-32)*Jnfl OOS4 t ]m 0 9 

where ♦ sinY 0 cosa 0 «const Is an Integral which corresponds to 
cyclic coordinate 0. We consider that Af*l. System (5) permits the 
following solutions: 
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XQ.*/*-m* V i cm*. . v («) 

c^r.,0, »-Q^*Ji'[«*-U)*« t jr.)en..a i XI) 

Go&y $ *0 9 s*f7^ • 0 (g) 

Angle v between the Ox^ axle of the natural rotation of the satel- /8 
llte and the current radius vector of the satellite center of mass 
relative to the center of mass of the earth Is determined by the re- 
lationship 

cog i> •eosfi MV *, , 

the angular velocity of the natural rotation Is determined by the 
equation 

jB»J20*s//f/ 0 eauc 0 . 

By using Integral (4) and the motions of the satellite linearised In 
the vicinity of the stationary solutions of the equations, both suf- 
ficient and necessary conditions of stability of these solutions can 
be obtained [4]. 

Solution (6) 

The sufficient condition of stability Is satisfied If. 

A-KO. 

The necessary conditions of stability are satisfied In the following 
two regions: 

J- 1* t\ 

The condition of existence of solution (6) has the form 

Uo 0 kH-3A| 
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The sufficient and necessary conditions of stability of the 
solutions oolnoldo and hava tba fora 

X-1>0, slno 0 P0. 

To thasa conditions must ba addad tha condition of existence of so* 
lutlon (7) 

I 


Solution (8) 


Let Yq*V2. 
tha form 


Tha sufficient conditions of Its stability hava 
XQ 0 -1>0, XQ 0 +3X-H>0, 


and tha naoassary conditions of stability arc /9 

(X& 9 -i)(XQsaM)*C, 


M. Satellite Motion Close to Stationary Rotation (6) 

We will investigate tha foread 2* periodic solution of system 
(2) which satisfies tha boundary conditions 


by solving boundary value problem (2), (9)* Va investigate boundary 
value problem (2), (9) by the Poincare small parameter method [53* 

Va use solution <$) as tha generating solution. For tha other gen- 
erating solutions examined in this work, such an Investigation Is 
oarrled out similarly. Va use vector notation to shorten the writing 
[93. Ve Introduce vector Z«(V^,V 2 »fi»Y»e)^* *nd we define function 
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P u t *f)£ £* so that system (2) and boundary oonditlona (9) 
could be written In the form 


and 


t(£l)-Kd“0 


(a*) 

( 9 ') 


roopoctlvoly. Ut t 

f(u,a,\ty,Z / be the solution of eyeteu (2* ) with the Initial 
conditions ?(0,a,l f kf)”as(a^, . . *-) . Boundary value p rob lea 

(2')» (9') can then be written 

g(o.i, kf)t iim i, (l0) 

We will consider relationship (10) as an equation relative to a. If 
kg«0, this equation permits the solution aiCa^. • • • » a^) 1 where 

4*4 4 — 4 v>, < n > 


Because of the analytical nature of the right side of system 
(2*) with respect to Z and kg, with sufficiently small kg and |Z-a|, 
function g(a,l,kg) analytically depends on kg, a In the vicinity of 
point kg»0 and a-a. If £ 12 . 

t * 0, ( 12 ) 

according to the theorem of the Implicit function, with sufficiently 
small kg, Eq. (10) has the unique solution a-2(kg,l), which depends 
analytically on kg and satisfies the condition 2(0,i)»a* In this 
case, boundary value problem (2')» (9*) has the unique solution 

f • i, kj ), (13) 

which depends analytically on kg in the vicinity of point kg-0 and 
coincides at this point with stationary rotation (6) with o Q -0. 
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We win investigate solution (13) in the fcrs of an Integral 
power series of parameter kg 

0-4*44,*-, (1*) 

> d*Gk*4d*u)+... . 

with 2t periodic coefficients with respect to u. The equations in 
the variations for stationary solution a have the fora 

A -XoetfisQ. ( 15 ) 

aA» 0, Af*AW r sinf 0 AJ> JU*dW s t$frjLd?. 

me characteristic equation of system (15) is separated into the 
equations p»0 and 

P i +op*+f*0 i ( 16 ) 

where am ?-£La-$2/P3)sbfy, f Por determlnatlon of 

the forced solution of system (2), we substitute series (14) in sys- 
tem (2) , and we equate the terms with the same powers of kg. We 
obtain a series of systems, the general form of which is the following 

%**(*-%**&«* '3a;*, . (17) 

'it* WicSiri'it.f * T$*»r»fK+&s t g. 


Here, S 1 k (i»l, ... 5; k«l,2, . . .) are some functions W^ 0 , .... /II 

W 1 k-1* • • *i °o* • • •» a k-l 5 u * S l,o“°* * 7 * t9m (17) » t0 wlthln 

the designation of the variables with k*0, coincides with system (15). 

In the solution of system (17), we will use the results of [63* Let 
the solution of system (17) be found up to k-1 inclusive. This solu- 
tion depends on integration constant 0 k _^. The equation for Q k is 
separated out, and it can be integrated. Then, 
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(t)e/r*Ji' M . 

The condition of 2tp riodicity of the function 

?S v (t)enr.0 ( 18 ) 

permits determination of the value of h ^ . In particular, to find 

function Q^(u), we have the equation 


According to Eq. (18), function fi^(u) is 2ir periodic If sinY 0 *0 or 


7*3^) 



To within 0(kg), the first condition corresponds to the motion of the 
satellite in the orbital plane. From the second condition, we find 
the region of existence of generating solution (6) in the plane 
(sin 2 i, X). The boundaries of this region are determined by curves 
X^(i) , X 2 (i), which are assigned by the expressions 




V-cat* 


J)+5sir?i 






The equality cos 2 y 0 "0 (cos 2 y 0 *1) 1» satisfied along the curve 
X^i) (X 2 (i)). The region within which the inequality 0<cos 2 y 0 <1 
is satisfied is crosshatched in Fig* 2. Curves X^(l) and X 2 (1) 
intersect at point P with coordinates 
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Subsequently In this section* we consider that y Q and consequently /12 
Qq are determined by Eq. (19)* 

Expansion of the following functions In a Fourier series Is 
subsequently required: 

Stnu _ 2 fp , 

'VTtoiVZ&Tih " sm(*n.,) u , 


*$ m 3*>* i 






H* 3sin A i mf 


9e*"ce& 2m*. 


Here Is Introduced the notation 


i s/nd 


where 


The expression for can now be written. It has the form 




nyt^sx/Ji 


The value of 0^ will be determined below. By substituting the ex- 
pression found for 8^(u) In the equations of system (17) for 
W 2 y^* a v we obtain a system of linear heterogeneous equations 
with periodic free terms. The corresponding uniform equations coin- 
cide with the equations of system (15)* In which the equation Afl«0 
should be excluded* Aft*0 should be set, and the variables should be 
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ignated. We assume that Eq. (16) does not have a root of the 
p»k«£l with any whole k, i.e. , # 

**-K*a+ 6? 0, (20) 

eterogeneous system In question then has the unique 2 * periodic /13 
i° n W 1#1 (u), W 2 i^)» a^(u). This solution can be found 

|e form of trigonometric series of the form 

I _ 

£ 0 {t*a*nu+£ n 's;»nu) (21) 

jso can be shown that, upon satisfaction of condition (20), such 
s exist and are unique. Prom the condition of 2 * periodicity of 

Ion n 2 (u), we find An example of construction of the 

ion in the form of trigonometric series in explicit form will be 
iited in Section 6. The solution for arbitrary k is constructed 
similar manner, if all solutions to some k-1 Inclusive are found. 

2 

We find the values of X* A (sin 1) at which condition (12) Is 
ted. It can be shown that J«0 when and only when system (15) 

>s a nontrivial solution which satisfies the boundary conditions 

by writing out the general solution of system (15), we find that 
sen and only when solution (16) has the root p«k,£T with some 
k. Thus, condition (12) is equivalent to condition (20). Cal- 
jLons have shown that condition (20) is violated at k*0 on the 
curve (Pig. 2). At k*l, condition (20) is violated on the res- 
s curves designated by the dashed lines. If k>l, the resonance 
s pass through point p, but they lie outside the region of exlst- 
Df generating solution a. 

Leal Study of Satellite Motion Generated from Solution (6) 

for arbitrary values of parameter kg, we construct solution (13) 

13 
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numerically, by solving boundary value problem (2), (9), and we in- 
vestigate its dependence on parameters kg, X. The solution of this - 
boundary value problem is reduced to solution of systom (10). Here 
and subsequently, system (10) is considered a system of equations which 
defines the curve in space !r(a, kg) with X*oonst [7]. System (10) 
was solved numerically by the methodtof Newton, in which both a and 
kg were refined at each step. For calculation of functions 

g(a,:kg), *8) , which are used in the method of Newton, system 

o 

(2) and the system in -variations corresponding to it were Integrated 
in the Interval 0«u<2ir. The solutions of boundary value problem (2), 

(9) presented in Fig. 5 were found by this method. Figure 5 also /It 

presents the dependence of a^, . . . , a^ on kg with X- const, l»*/2. 

It is easy to obtain the other curves by means of substitution (3)« 

Here and subsequently, the number beside the curve in the figures 
designates the value of fixed parameter X. With i»w/2, the right 
sides of the equations in system (2) are v periodic functions of the 
true anomaly , which permitted the substitution 2rMr to be performed 
in boundary conditions (9) and permitted restriction to the integra- 
tion Interval 0«u{jr* 

For convenience, we will call this method of construction of 

g 

curves in space R (a, kg) extension by parameter kg. Three types of 
solution of stationary rotation (6), obtained by extension by kg, 
can be distinguished. In the Interval V7<X<2/3, solution (13) is 
extended right up to merger with the plane solution which describes 
the motion of the axis of symmetry of the satellite in the orbital 
plane. Branching of the solution within the region of its existence 
does not occur. In the Interval 2/3<X<8/ll, solution (13) also is 
extended right up to merger with the plane solution, but branching of 
this solution within the region of its existence occurs. In the in- 
terval 8/ll<X<V5, merger of solution (13) with the plane solution 
does not oocur. This solution "escapes” to the region of larger 
values of kg. Interval V7<X<4/5 of existence of * periodic solutions 
is broken down into the Intervals Indicated by points X»2/3 and X*8/ll, 
from which the resonance curves for k*l originate. 
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Amplitude curves of the solutions obtained art presented In 
5» We understand amplitude here to be the quantity 

( 22 ) 




where e , r are unit vectors along the Oz axis and the axis of sta- 
tionary rotation of the satellite respectively. The position of the 
latter In space Is determined by relationships (19) and «q«0. 


The stability of the solutions obtained was Investigated In the 
following manner. The system In variations along solution (13) which 
corresponds to system (2) was integrated. Roots pj, . . . , of the 
characteristic equation for the system In variations were calculated. 
Degree of stability X„ of the resulting periodic solution, which de- /15 
termlnes the response speed of the system, was calculated by the 
equation 

J, enmoxljsj 

The condition X >0(X_*0) corresponds to a stable (unstable) solution. 

8 8 

Sections of the X_ curves which X >0 are presented In Fig. 6 with 
various X. The nature of the roots which determine X g changes at the 
break points of the curves. The curves marked with hachures In Fig. 

5 Csic] correspond to stable solutions. 


We extend the periodic motions of the satellite plotted with 
i*ir/2 by parameter 1 for kg*kg f >0. It can be proved by the Poincare 
small parameter method that, because of analytical nature of the right 
side of system (2) with respect to Z and 1, with sufficiently small 

//- /, /§.%/<* / t function g(a.,i,kg) depends analytically on 1, a 

in the vicinity of point i»ir/2, a»i Q , kg«kg*. Here, Z Q is solution 
(13) with i*ir/2, u*0. If 

J •a/til bO, 


according to the Implicit function theorem, with sufficiently small 
1 1- 1 / 2 1 , Eq. (10) has the unique solution 
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<),<,%), ( 23 ) 

which depends analytically on 1 In the vicinity of point i»V2 and 
coincides at this point with solution (13) with i-w/2, which was 
constructed above. This sane solution can be obtained by extension 
of solution (13) by parameter kg with !■!*. This method of extension 
of periodic solutions in the (i, kg) plane was described in detail In 
C83. 


Solution (23) was constructed numerically. The results of the 
calculations, which were performed for kg«0.2 and several values of 
X, are presented in Fig. 7. Curves of the initial values of phase 
variables 6 m and X 8 are represented here by the solid lines. The 
curves of the tabulated values for the solutions which are not char- 
acterlzed by rotation, but by oscillations of the satellite around 
its axis of symmetry, are represented by the dashed lines. Such 
solutions were constructed in [1]. Solution (23) exists right up to 
the point of merger with the solution, the curves of which are desig- 
nated by the dashed lines. 

The explicit form of periodic motion at X*0.63» kg*0.2, i»1.37 
is presented in Fig. 8. Curves of the phase variables, angle 6 and the 
trace of the Oz axis on a unit sphere which surrounds point 0 are 
presented here for 0<nf2, n is the number of orbits. Angle 6 is de- 
termined by the expression (see Eq. (22)) 

z). 

Where possible, the curves of the corresponding stationary solution /16 
a are designated by dashed lines and the symbol (•). The arrow on 
the curve in the (a,y) plane indicates the direction of motion of the 
trace of the Oz axis with Increase in u. The points on the curve are 
0.1 orbit apart. The values of y Q and were determined from Eq. 

(19) and (6). 
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By using the algorithm and notation of Section 4, we will seek 
a solution of system (2) in the form of integer power series (14) of 
parameter kg. As solution a of system (10) with kg*0, we seleot sta- 
tionary rotation (7): 

<£-*. (24) 

The equations in variations for the stationary solution selected have 
the form 

J iW,*Xtirrt.*0, (25) 

cos*!, > tiL* 


The characteristic equation of system of equations (25) Is broken down 
into the equation p-0 and 

. 0 . (26) 

For determination of the forced solution of system (2), we substitute 
series (14), where the first terms are determined by Eq. (24), in Eq. 
(2), and we equate the terms for identical powers of kg. We obtain a 
series of systems, the general form of which is the following 

* \* *3(1- , 

V Vi,* * - Wyt +Jt&t/7* 0 • Q. jc , 

Here, T t k (i-1, . . . , 5) are some functions W x 0 , . . ., W x k-1 , 
... , a*, . . ., a^, u, T t o »0. Similarly to the way it was done 
in Section 4, for determination of the constant value of a 0 and con- 
sequently, the value of cosa 0 , we write the equation for fi 1 (u) 

^ co $*<*(6ff Sin* e (s2h t1 hg i y/*<'et*( 0 
*■ h u *i4)jhg$)] • 
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The condition of 2w periodicity of function Q^(u) 

fT Am( F)srm0 

is reduced to the equalities coso 0 *0 or 


01 



( 28 ) 


The equality coscXq^O Is satisfied with any permissible X and 1. It 
corresponds to orientation of the axis of symmetry of the satellite 
along the velocity vector of its center of mass. Investigation of 
this generating solution can be carried out within the framework of 
another sequence of rotations a, 6, y» that presented in [1] for 
exanple. With X>1, such a solution is stable. From the condition 
Ogcosa 2 .$l, we find the region of existence of generating solution 
a (24) determined by equality (28) in the (sin 2 i, X) plane. The 
boundaries of the region are fixed by the X^(i), X 2 (l) curves, which 
are determined by the expressions 


7/;,, 


tV,. tfcss 


The equality cos 2 a Q »0 (cos 2 a 0 -l) is fulfilled along the X^d) (X 2 (i)) 
curve. The region bounded by the X^i), X 2 (i) curves is crosshatched 
in Fig. 3. The X x (i) and X 2 (i) curves intersect at point p, with co- 
ordinates * 0 , 6 #", 




JU 





By using the expansion of the functions Included in the right 
side of Eq. (27) in a Fourier series, we write the solution of Eq. 
(27) in the form — j*,, 

where 
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For determination of functions 1 , W 2 i* *i» ®i* we obtain a sys- 
tem of llnesr heterogeneous equations with periodic free terms. The 
corresponding uniform equations coincide with Eq. (25)* where the 
equation AQ*0 should be excluded and Afl"0 should be placed In the 
remaining equations. If Eq. (26) does not have the root p«kv^I at 
any Integer k, l.e.* 

k *0, ( 29 ) 

the resulting system of heterogeneous equations then has the unique /18 

2* periodic solution ^u), W 2 ^(u), y 1 (u)* o^u). It can be found 

in the form of a trigonometric series of the type of Eq. (21). A 
series for arbitrary k is constructed in a similar manner, if all 
solutions to some k-1 inclusive are found. 

The calculations showed that condition (29) is violated at k«0 
(k«l) on the (X » (1) ) curve. If k>l, the resonance curves pass 

through point p and lie outside the region of existence of generating 
solution a (24). 

For arbitrary values of parameter kg, construction of 2* periodic 
motions of the satellite is reduced to numerical solution of system 
(10) and extension by kg of the solution constructed In form (14), 
which Is close to stationary rotation (24). In the case of a polar 
orbit, the results of extension of solution (13) by kg. In the form 
of curves of functions a^, . . . , a^ with X*const, are presented In 
Fig. 9* The initial conditions of the solutions obtained from the 
solutions found by transformation (3) can be plotted by symmetrical 
representation of the curves (Fig. 9) relative to the corresponding 
axes. 


Investigation of the stability showed that all the solutions con- 
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struoted art unstable. The amplitude characteristics arc presented 
In Pig* 10. They wart datarminad by Eq. (22), whara tha position r 
ia datarminad by ra 1st ionahipa (28) and t 0 -t/2. 

6* Satellite Motion Cloaa to Stationary Rotation (8) 

By using tha algorithm and notation of Saetion A, wa soak solu- 
tion (13) of system (2), ganaratad from stationary solution (8), in 
tha form of aeries (14). Va select stationary rotation (8) as solu- 
tion a of system (10) with kg»0: 

*s m a (30 

The value of Q Q is subject to determination. The aquations in varia- 
tions for tha stationary solution selected have tha form 

AH', * -/f+J& 0 )AW M 

j (/tZSljjA Wf , dSl* % 

Aj'mAWf-A'l, J 


The characteristic aquation of system (31) is broken down into tha 
aquations p*0 and 

c'*p , [iXa..+tfT3}-z]*(lQs<t-3mSl.*<).C> ( 32 ) 

For determination of tha coefficients of series (14), we obtain a /19 
series of systems, the general form of which is easily described from 
Eq. (17), (31). The condition of existence of a 2w periodic solution 
of the equation 

A/u) . 03) 


is reduced to the expression 
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solution of Eq. (33) csn bt writtsn In ths form 

i, f£»*~**f 

If Eq. (32) doss not have ths root p«k/3i at any intsgsr k, i.e., 

k-m.i ( 3 " } 

' ^^nsn ths linsar systsm rslativs to variables k , W 2 lc> a k , with 
: ''^s=sTiodic frss terms, has a unique 2* periodic solution. Por k»l, vs 
v ~ - ~--ite this solution in ths following way 



6 


Wy *£[C ’(2n+t)Jn]*f '•(**+ 'J** 
^4.£fC n (Zn*t)-at n ]co*(2f»*)“t 

090 

af„ . i„(ll2»* f* fa* W&t *)*&&.* ')} , 







Prom the condition of existence of 2ir periodic function q 2 (u), 
"^nne equation for determination of which is not presented here be- 
— aause of the cumbersome form, we obtain n 1 °-0. The solution can be 
-constructed in the same manner, in the form of series for arbitrary 
if all solutions up to k-1 inclusive are known. 


Thus, in ths two prsosding ssctions, invsstigation with arbitrary j 
kg of solution (13) gsnsratsd from Eq. (30) was rsduesd to numerical 
-solution of systsm (10) with X»eonst. The calculation results are pre- 
sented in Fig. 11, in ths form of a^ curves for several values of X am 
' — l«w/2. Condition (34) is violated in the curves presented in Pig. 12. 
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